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ABSTRACT 

This paper reports the star count predictions of warped and flaring models 
of the outer Milky Way disk. These have been suggested as possible locations of 
the lenses responsible for the microlensing events towards the Large Magellanic 
Cloud (LMC). Three such models are investigated in detail and the theoretical 
predictions are confronted with Hubble Space Telescope (HST) star count data in 
7 low latitude fields (30° < |6| < 40°). If the warped and flaring disk population 
has the same characteristics as disk stars in the solar neighbourhood, then the 
models can be unambiguously ruled out. Metallicity gradients are well-known 
in disk galaxies and may cause the outer disk population to differ in colors 
and luminosity from that locally. This effect is studied using a simple ansatz 
for the variation in the color- magnitude relation with position and while it is 
shown to lead to better agreement with the star counts, upper limits on the 
contribution of the warped and flaring disk to the optical depth are still below 
the value measured towards the LMC. Only if the warp is very asymmetric or 
if the luminosity function changes strongly with Galactocentric radius can the 
models be made consistent with the starcounts. 

Subject headings: Galaxy: kinematics and dynamics - gravitational lensing - 
dark matter 
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1. INTRODUCTION 

The difficulty of explaining the origin of the microlensing events towards the Large 
Magellanic Cloud (LMC) with conventional models of the Milky Way has led to the 
suggestion that the outer parts of the stellar disk may be warped and flared (Evans, Gyuk, 
Turner & Binney 1998). Direct observational evidence on the extent of warping and flaring 
in the outer Milky Way is hard to gain, but models in which the stars warp and flare as 
strongly as the neutral gas can provide microlensing optical depths of ~ 1 x 10~ 7 . The 
aim of this paper is to examine this suggestion with a critical eye in the light of the new 
deep star count data that is available from the Wide Field Camera on the Hubble Space 
Telescope (e.g., Gould, Bahcall & Flynn 1998). The warp of the stellar disk is difficult to 
test directly. This is because searching closer than about 20° to the Galactic plane is difficult 
as the confusion due to regular foreground sources becomes too high. By comparison, the 
maximum warp angle is only about 8°, and then only along the ridge-line of the warp. An 
easier target is to search directly for the evidence of the flare that is assumed on top of the 
warp. It is the flare that is mainly responsible for increasing the microlensing optical depth 
towards the LMC to values broadly consistent with the observations. The flare must put an 
enhanced number of stars along the line of sight to the LMC (£ = 280°, b = —33°). In this 
paper, our aim is to confront the theoretical predictions of warped and flaring disk models 
with the star count data on low latitude fields with 30° < \b\ < 40°. 

Section 2 introduces the dataset under scrutiny. In Section 3, the simplest assumption 
is made - namely that the stellar population in the warped and flaring disk has the same 
luminosity function as that determined for the local disk. Section 4 considers more elaborate 
models in which the effects of metallicity gradients are taken into account, albeit crudely. 
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2. THE STAR COUNT DATASET 

The dataset is part of an on-going star count project being conducted with the Wide 
Field Camera on the HST (see e.g., Gould, Bahcall & Flynn 1996, 1997, 1998). The seven 
fields we have selected for scrutiny are listed in Table la. They have been chosen so that 
30° < |6| < 40° and so are optimum for our purposes of tracing the stellar flare at the 
lattitude of the Magellanic Clouds. The table also lists the I band saturation magnitude 
7 m i n and the detection threshold I max . The saturation magnitude varies from field to field 
as it depends on the exposure time of the shortest exposure used in producing the cosmic 
ray cleaned images. The detection threshold also varies from field to field - it is the faintest 
magnitude at which stars and galaxies can still be unambiguously distinguished. Full details 
of the methods of determining these limits are given in Gould, Bahcall & Flynn (1996). To 
be conservative, 1.5 magnitudes are added to the nominal bright limit and 0.3 magnitudes 
are subtracted from the nominal faint limit. On the faint side, this guarantees absolute 
safety with regard to star/galaxy separation. At the bright end, the 1.5 magnitudes are due 
to technical reasons concerning the pixel size. Only data points between these new limits 
are actually used in constructing the star count statistics. 

Table lb gives the star count data in the color range 1.52 < V — I < 2.26 in the seven 
fields. This selection function picks out M dwarfs, as discussed in Gould, Bahcall & Flynn 
(1996). The M dwarfs should be representative of the disk, and so should presumably 
participate in any outer flaring or warping. Table lb shows that each field contains only a 
few stars. The numbers of stars in each of seven magnitude bins is recorded. The bins are 
of 1 magnitude extent. The / band magnitude centroids of the bins are indicated at the top 
of each column. 
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3. WARPED AND FLARING DISK MODELS 

The density of the warped and flaring disk model as a function of galactocentric radius 
R and height z above the Galactic plane is: 

p(R,z) = p exp(-R/R d )\(l - /) exp(-|z - z w \/z mck ) + /sech 2 (|z - z w \/z thin )\. (1) 

The overall normalisation p is chosen to reproduce the column density of the thin and thick 
disks at the solar radius (R Q = 8kpc) of 46±3M pc -2 , in agreement with recent dynamical 
estimates (Gould 1990, Holmberg & Flynn 1998). More important from the standpoint of 
comparison with star counts is the normalisation of the disk with the contribution of the 
gas (~ 13M pc~ 2 ) and the disk white dwarfs (~ 4M pc~ 2 ) excluded. This column density 
is taken as ~ 3OM pc~ 2 . (Let us remark that this value inferred from dynamics is perhaps 
some ~ 15% higher than suggested by Gould, Bahcall & Flynn (1997) on the basis of 
starcounts. These authors reckon that the total column density in stars is greater than the 
column density in M dwarfs by a factor of 2.1, leading to a value of ~ 26M pc~ 2 ). 

The vertical profile of the thin disk follows roughly a hyperbolic-secant squared profile 
whereas the thick disk is modelled by an exponential profile in height above or below the 
midplane of the warp \z — z w \. The fractional contributions of the thin and thick disk to 
the midplane density are denoted by / and 1 — /. The sun lies on the line of nodes of the 
warp, which is described by the equation: 

z w (R,9) = 0.214 x (R-R e ) sin9 (2) 

and sets in for radii R greater than the solar position R Q . Here, 9 is the longitude measured 
about the Galactic Center from the line of nodes. Observationally, the stellar warp cannot 
be detected to large radii. Our equation is suggested by the behaviour of the gas warp with 
one important difference. The gas warp bends back to the Galactic plane in the southern 
hemisphere for R > 15kpc (Burton 1988), but the amplitude of the stellar warp is assumed 
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to continue increasing linearly. The scaleheights of the thin and thick disks flares by a 
factor F between the solar position and 25 kpc and so are described by; 



^thin — ^thin,0 



R-R P 



( F ~ — TT~ , ^thick = ^thick,© 1 + (F - 1) " . (3) 
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For comparison, the scaleheight of the neutral gas is known to flare by a factor F of roughly 
ten over this distance (e.g., Wouterloot, Brand, Burton & Kwee 1990). Table 2 gives the 
scale length, flaring factor, relative normalization of the thin and thick disks and the scale 
heights for the 3 models representative models we consider here. For all models the disk 
is truncated at 25 kpc. Also shown is the optical depth towards the LMC. Effects such as 
ellipticity of the disk, or lopsidedness, which can increase the optical depth further are not 
considered. 

Given the density law, a luminosity function $ is needed to make star count predictions. 
The M dwarf luminosity function locally is recorded in Fig. 2 of Gould, Bahcall & Flynn 
(1998). Between 8 < M v < 10.5, the V band luminosity function can be approximated by 
the linear fit 



$ = 0.0633 x 



3.16(M„ - 8.) + 1.2 , (4) 



where $ is normalised to unit density. As the star counts are known in the / band, but the 
luminosity function is readily available in the V band, we will need to convert using the 
disk color-magnitude relation (Gould, Bahcall & Flynn 1998) 

M v = 2.89 + 3.37(V - I) (5) 

and assuming vanishing extinction so that V — I — M v — Mj. The results for our three 
models are given in Table 3. It is clear that if the characteristics of the outer Milky Way 
disk are the same as those of the local solar neighborhood, all three models are strongly 
ruled out. This conclusion has been confirmed by Holmberg (1998, private communication) 
in a starcount analysis of a somewhat similar model - namely the modified Bahcall-Soneira 
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model discussed in Holmberg, Flynn & Lindegren (1997) on which a flare and warp have 
been superposed. 

Neglecting the contribution of the spheroid, we perform a maximum likelihood fit of 
the overall normalisation of the stellar disk to the star count data. The fit is best - though 
still poor - for a disk surface density (participating in the flare) of ~ 2.5M pc~ 2 . This gives 
an upper limit to the optical depth available from the warped and flaring disk of 3 x 1CT 9 . 

4. THE EFFECT OF COLOR GRADIENTS 

The assumption that the disk stellar population is constant with radius is attractive 
from the point of view of simplicity. Unfortunately it is hardly likely to be correct. Color, 
metallicity and even mass function gradients are well established in disk galaxies (Pagel 
1998, Taylor 1998). Spheroid stars in this color range are typically 2.5 magnitudes fainter 
than disk stars of the same color (Monet et al. 1992). All this suggests investigation of 
color-magnitude relations whose properties change with Galactocentric radius such as 

M v = 2.89 + 3.37(V - I) + l.ffi\(R - R e ) (6) 

Here, A is the metallicity gradient in dex kpc -1 . We have assumed that the change in 
magnitude between disk and spheroid is AM V = 2.5 and the corresponding change in 
metallicity is AZ = 1.5. Following the findings of Hippel et al. (1996), we assume that the 
luminosity function for the range we are interested in is independent of metallicity. 

The results of adding a metallicity gradient can be seen in Table 4 which shows the star 
counts expected for model A with A = 0.15 dex kpc~ x . Clearly, the number of stars expected 
has been substantially reduced. However, performing the same maximum likelihood fit as 
previously, we find that the upper limit on the optical depth has only improved ~ 1CT 8 , still 
far below the measured optical depth towards the LMC This is perhaps hardly surprising. 
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It is well-known that the deep star count data are consistent with a spheroid whose density 
falls off like r -3 - 5 and which contributes at most a few percent to the optical depth. If a 
warped and flaring disk is to be consistent with the same star count data, its optical depth 
must be nearly the same to within a geometrical factor of a few. 

5. CONCLUSIONS 

The warped and flaring disk models of the outer parts of the Milky Way cannot provide 
a substantial contribution to the microlensing optical depth and still be consistent with the 
starcount data. This statement needs two qualifications. First, the starcounts presented 
here are not exactly in the direction of the Large Magellanic Cloud (LMC). Strictly 
speaking, it is therefore the assumed azimuthal variation of the warp and flare that is being 
directly tested. It remains possible that a strongly asymmetric warp (perhaps generated 
directly by the LMC) can be made consistent with the starcounts. Second, the luminosity 
function or mass function may vary strongly with radius and this may allow the lenses to 
remain dark enough to evade detection even in very deep fields. This idea has sometimes 
been regarded as unappealing when measured against the yardstick of theoretical simplicity, 
though it has recently been resuscitated by Kerins & Evans (1998) and Binney (1998). In 
particular, there may be two modes of star formation controlled by metallicity. It might 
be that at low metallicity, only low mass objects are formed, whereas at higher metallicity, 
mainly viable stars are produced. 

Many researchers have suggested that the microlensing events seen towards the LMC 
may occur not in the Galactic halo but in intervening structures. Possibilities include 
not just the warped and flaring disk of the Galaxy, but also stellar conglomerations in 
front of the Magellanic Clouds (Zaritsky & Lin 1997, Zhao 1998) and hitherto undetected 
populations in the outer parts of the Galaxy (Gates, Gyuk, Holder & Turner 1998). If 
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these populations are traced by at least some bright stars, then they should be visible in 
deep starcounts towards the LMC. So, all these suggestions are very vulnerable to direct 
observational testing both with the Hubble Space Telescope and the new generation of large 
telescopes (such as the Very Large Telescope). 

GG thanks the Theoretical Physics Department for hospitality during a working 
visit to Oxford University, the Department of Energy for partial support under grant 
DEFG0390ER40546 and Research Corporation. NWE is supported by the Royal Society. 
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to use their star count data in advance of publication. We also thank Johan Holmberg for 
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Field No 


i 


b 


-^min 


-^max 


4 


179.828° 


-32.146° 


17.2 


23.9 


36 


56.716° 


34.251° 


18.3 


24.2 


38 


56.726° 


34.248° 


18.3 


24.2 


71 


149.784° 


34.678° 


18.0 


23.8 


72 


149.784° 


34.712° 


18.4 


24.0 


94 


303.348° 


33.66° 


18.3 


24.1 


140 


77.018° 


35.78° 


18.1 


24.2 



Table la: This table lists the seven fields under scrutiny. The first column gives the 
(arbitrary) field identifier and the next columns give the location in terms of Galactic 
longitude and latitude. The final two columns give the limiting I band magnitudes. 
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Table lb: This table lists the number of stars with colors satisfying 1.52 < V — I < 2.26 in 
the seven fields. There are seven histogram bins each of 1 magnitude extent. The I band 
magnitude centroids of the bins are indicated at the top of each column. The starcounts in 
the final column may be incomplete because of the upper limit on the I band magnitude. 



Model R d F f z thhl:Q z thick:Q r 



A 3.5 10 0.75 350 700 3.8 x 10~ 8 
B 4.5 10 0.90 350 1000 4.4 x 10~ 8 
C 4.5 5 0.75 350 700 2.3 x lO" 8 

Table 2: This table gives the parameters for the three warped and flaring disk models 
examined in detail. They are the scalelength Rd and the scaleheights of the thin disk z t hm,e 
and the thick disk z t hi C k,© at the solar radius. The flare parameter F describes the factor 
by which the scaleheights increase at a Galactocentric radii 25 kpc. The final parameter / 
describes the relative proportion of the thin disk to the thick disk. For each model, the 
microlensing optical depth r is given in the final column. 
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Field No Model 18 19 20 21 22 23 24 
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Table 3: This table lists the expected number of stars with colors satisfying 1.52 < V — I < 
2.26 in the seven fields for each of the three warped and flaring disk models. 
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Field No 18 19 20 21 22 23 24 
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Table 4: This table lists the expected number of stars with colors satisfying 1.52 < V — I < 
2.26 in the seven fields for model A assuming a metallicity gradient of 0.15 dex/kpc from 
the solar circle. 



